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Abstract: The difference between vacuum energy of quantum fields in Minkowski space 
and in Friedmann-Robterson- Walker universe might be related to the observed dark energy. 
The vacuum energy of the Veneziano ghost field introduced to solve the ^7(1)^4 problem in 
QCD is of the form, H + 0(H 2 ). Based on this, we study the dynamical evolution of a 
phenomenological dark energy model whose energy density is of the form aH + /3H 2 . In 
this model, the universe approaches to a de Sitter phase at late times. We fit the model 
with current observational data including Snla, BAO, CMB, BBN, Hubble parameter and 
growth rate of matter perturbation. It shows that the universe begins to accelerate at 
redshift z ~ 0.75 and this model is consistent with current data. In particular, this model 
fits the data of growth factor well as the ACDM model. 
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1. Introduction 

The accelerating expansion is still a mystery of modern cosmology since its discovery in 
1998 [1], and a new energy component called dark energy (DE) is needed to explain this 
acceleration expansion within the framework of general relativity. The simplest model of 
DE is the cosmological constant, which is consistent with all observational data, but it 
faces with the fine tuning problem [2]. Instead, many alternative DE models have also 
been proposed [3, 4, 5, 6, 7, 8, 9], but almost all of them explain the acceleration expansion 
either by introducing new degree(s) of freedom or by modifying general relativity. 

Recently the so-called QCD ghost dark energy has been proposed in [10, 11, 12]. The 
Veneziano ghost field plays a crucial role in the resolution of the U(1)a problem in QCD 
[13]. The ghost field has no contribution to the vacuum energy density in Minkowski space- 
time, but in a curved spacetime, it gives rise to a vacuum energy density proportional to 
A 3 QCD H [14, 15, 16, 17], where A QC d is QCD mass scale and H is Hubble parameter. Note 
that in this ghost dark energy model, there are no unwanted features such as violation of 
gauge invariance, unitarity, causality etc. [10, 11, 14, 15, 16, 17]. In fact, the description 
in terms of the Veneziano ghost is just a matter of convenience to describe very compli- 
cated infrared dynamics of strongly coupled QCD. The Veneziano ghost is not a physical 
propagating degree of freedom, one can describe the same dynamics using some other ap- 
proaches (e.g. direct lattice simulations) without using the ghost. Therefore the Veneziano 
ghost field is quite different from those ghost fields in some dark energy models in the 
literature, there those ghost fields are real physical degrees of freedom, introduced in order 
to have the equation of state of dark energy to cross —1. On the other hand, the vacuum 
energy is generally expected exponentially suppressed because QCD is a theory with a 
mass gap. However, this issue is elaborated in some details in [14, 15, 16, 17], there it has 
been convincingly argued that the complicated topological structure of strongly coupled 
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QCD will lead to a vacuum energy density with an inverse linear scale of manifold size in 
a nontrivial background. The power law behavior is also supported by recent lattice result 
[15, 18]. Very recently it has been shown that this behavior is also got supported from the 
holographic description of gauge field [19] . 

Because this model is totally embedded in standard model and general relativity, one 
needs not to introduce any new degrees of freedom or to modify Einstein's general relativity. 
In this model, the energy density of DE is roughly of order Aq CD H, with Aqcd ~ lOOMeU 
and H ~ 10~ 33 eU, so A 3 QCD H gives the right order of observed DE energy density. This 
numerical coincidence is impressive and also means that this model gets rid of fine tuning 
problem [10, 11]. The model parameters have been fitted recently by observational data 
including Snla, BAO, CMB, BBN and Hubble parameter data [12]. It shows that this 
model is consistent with those observational data. 

On the other hand, it is convincingly argued that the contribution of zero-point fluctu- 
ations of quantum field to the total energy density should be computed by subtracting the 
Minkowski space result from the computation in a FRW space-time [20, 21]. Usually the 
difference, H 2 A 2 , between the vacuum energies in Minkowski space and in the FRW space- 
time, is absorbed into a renormalization of Newton's constant G, here H is the Hubble 
constant of the FRW universe and A c is the cutoff. However, it is true only under the as- 
sumption that the vacuum expectation value of the energy-momentum tensor is conserved 
in isolation [21]. In that reference the authors investigated the role of this term as early 
dark energy in the evolution of the universe. 

Notice the fact that the vacuum energy from the Veneziano ghost field in QCD is of 
the form H + 0(H 2 ), see for example, [19], while in the previous works on the QCD ghost 
dark energy model, the leading term H has been considered only. Having considering the 
study in [21], one may expect that the subleading term H 2 in the ghost dark energy model 
might play a crucial role in the early evolution of the universe, acting as the early dark 
energy. 

Based on the QCD ghost dark energy, in this work we therefore investigate a phe- 
nomenological model with energy density pde = otH + f3H 2 . For other motivations to 
consider this form see [22], there it is argued that this form of varying cosmological con- 
stant could be a possible candidate to solve the two fundamental cosmological puzzles. We 
study the cosmological evolution of this DE model, fit this model with current observa- 
tional data and give constraints on the model parameters. Besides, it is worth noticing 
that a varying DE should have some effect on the evolution of the matter perturbation, 
so we study the first order perturbation to the matter density and fit this model with the 
data of linear growth factor. 

The paper is organized as follows. In Section 2 we study the dynamical evolution of 
the DE model. In Section 3, we fit this model with current observational data and discuss 
the fitting results. The data used are Union II Snla sample [23], BAO data from SDSS 
DR7 [24], CMB data (R, l a , z*) from WMAP7 [25], 12 Hubble evolution data [26, 27] and 
Big Bang Nucleosynthesis (BBN) [28, 29], and we also study the effect of the DE on the 
linear perturbation of matter density. We summarize our work and give some discussions 
in Section 4. 
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2. Dynamics of QCD Ghost Dark Energy 



To study the dynamics of the DE model, we consider a fiat FRW universe with three energy 
components: matter, DE and radiation. In this ghost DE model, the energy density of DE 
is given by pde = oiH + f3H 2 , where a is a constant with dimension [energy] 3 , roughly of 
order of Aqq D where Aqcd ~ lOOMeV is QCD mass scale, and /3 is another constant with 
dimension [energy] 2 . For convenience, we define 7 = 1 — and use this throughout the 
paper. 

Arming with this DE density, the Friedman equation reads 



H 2 = 8 ^-( a H + Pm + p r ), (2.1) 
37 

where p m is energy density of matter, whose continuity equation gives 

Pm + 3Hp m = =^ p m = PmOCl 3 . (2.2) 

and p r is energy density of the radiation, whose continuity equation gives 

pr + 4Hp r = =>■ p r = PrO«~ 4 - (2.3) 

We have set ao = 1 and the subscript stands for the present value of some quantities. 
Solving the Friedman equation, we have 



4ttG UttG \ 2 8itG , 8irG . . lS 

H ± = —a±^ a j + ^P^oa- 3 + — p r0 a^. (2.4) 

There are two branches, H+ represents an expansion solution, while H- a contraction one. 
We neglect the latter since it goes against the observation, and for simplicity, write H + as 
H in what follows. 

Expressed with fraction energy density of matter and radiation, £l m o, ^ro, Equation 
2.4 gives an important constraint among these parameters: 

(7 - ^mO - ^ro)#0 = ~^~ a - ( 2 - 5 ) 

Further, Equation 2.4 can be rewritten as 



H (,)^o(K + ^ + ""°' 1 + ^ 7 +ii -°" + " 4 j. (2.6) 

where k = (1 — (O m o + Q-yo)/l)/2 an d z is the redshift, z = 1/a — 1. We can see from 
(2.6) that the universe approaches to a de Sitter phase with Hubble parameter 2kHq at 
late times, while it is dominated by matter and radiation terms at early times. 

Note that there exist many papers focusing on the coupling between the time-dependent 
vacuum energy and matter [30]. But we do not consider here such coupling in our anal- 
ysis. Namely in our discussion, the DE, matter and radiation are separately conserved. 
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In that case, the corresponding time evolution equation for the matter density contrast 
D = Spm/pm is given by: 

D + 2HD - 4irGp m D = 0, (2.7) 

where the over dot denotes the derivative with respect to the cosmic time. In terms of the 
growth factor [31], Equation 2.7 can be rewritten as 

-(1 + z)H{zff z + 2H{zff + H{zff - (1 + z)H(z)^f = + z) \ (2.8) 

where the growth factor / is defined as / = — (1 + z) ^ D ■ in general, there is no analytical 
solution to Equation 2.8, and we need to solve it numerically. But it is very interesting 
that the solution of the equation can be approximated as [32] 

/ = tt m (z)\ (2.9) 

and the growth index A can be obtained for some general models as 

where w is the equation of state of DE. For the case with 1 — Q m (z) being between zero and 
0.8, the accuracy is better than 1%. For the ACDM model, the approximation f(z = 0) = 
n^ + O A0 (l + Q m0 /2)/70 can be made [33]. But in our analysis, instead of parametrization 
of A [34], we will solve the Equation 2.8 numerically, by setting the initial condition f(z = 
0) = /o, where /o is a free parameter to be constrained by observational data. 

3. Data Fitting 
3.1 Model 

In order to fit the model with current observational data, we consider a more realistic 
model which includes DE, Cold Dark Matter, radiation and baryon in a flat FRW universe 
in this section. In this case, the dimensionless Hubble parameter can be written as, 



e = M. = k+ . / 2 + fU(i + ^ + o r »(i + ^ (31) 

Hq y 7 

where the energy density of baryon and Cold Dark Matter are always written together as 
SIdmo + ^60 = fimo, and Obmo, ^&o> ^rf) are present values of dimensionless energy density 
for Cold Dark Matter, baryon and radiation, respectively. The energy density of radiation 
is the sum of those of photons and relativistic neutrinos 



Q r0 = Oyo (1 + 0.2271iV n ) 

r 2 7 



where N n = 3.04 is the effective number of neutrino species and ft 7 o = 2.469 x 10^ 5 /T 



for T cmb = 2.725K (h = H /100Mpc ■ km ■ s^ 1 ). 

We will choose h, 7, tlbo and Q m Q (and also fo when we consider the growth factor) 
as free parameters of the model in the following data fitting. This relation Equation 2.5 
implies that there exists a strong degeneracy among h, 7 and f2 m g. 
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3.2 Observational Datasets 

We fit our model by employing some observational data including Snla, BAO, CMB, Hubble 
evolution data, BBN and the data of growth factor. 

The data for Snla are the 557 Uion II sample [23]. x 2 sn f° r Snla is obtained by 
comparing theoretical distance modulus nth( z ) = 51og 10 [(l + z) J dx/E(x)\ + no (no = 
42.384 — 5 log 10 h) with observed nob °f supernovae: 

Xsn = }_ J ^ • (3-2) 

To reduce the effect of no, we expand x 2 n with respect to no [35] : 

xL = A + 2Bno + Cnl (3.3) 



where 



-4 = £ 

i 



[nth(zi-,no = o) - nob(zi)]'' 

<J 2 (Zi) 

nth{zi-,no = o) - nob(zi) 



<T 2 {Zi) 



c 



y— 



(3.3) has a minimum as 

Xsn Xsn,min ^ ^ l^i (3-4) 

which is independent of no- in fact, it is equivalent to performing an uniform marginaliza- 
tion over no, the difference between x^n an d the marginalized x 2 n ^ s J us t a constant [35]. 
We will adopt x 2 n as the goodness of fit between theoretical model and Snla data. 

The second set of data is the Baryon Acoustic Oscillations (BAO) data from SDSS 
DR7 [24], the datapoints we use are 

, r s (z d ) 

and 

, _ r s{Zd) 
"0.35 



where r s (z ( i) is the comoving sound horizon at the baryon drag epoch [36], and 



D v (z) 



dx \ 2 z 



1/3 



H(x)J H{z) 

encodes the visual distortion of a spherical object due to the non Euclidianity of a FRW 
spacetime. The inverse covariance matrix of BAO is 

x ( 30124 -17227 \ 
M,bao ~ I _i7227 86977 J ' 
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The x 2 °f the BAO data is constructed as: 



Xbao 



V 1 "c^ 1 V 
1 M.bao ' 



(3.5) 



where 



y 



do.2- 0.1905 
do.35- 0.1097 



The third set of data we use are CMB datapoints (R, l a ,z*) from WMAP7 [25]. z* is 
the redshift of recombination [37], R is the scaled distance to recombination 

dz 



R — V^hnO 

and l a is the angular scale of the sound horizon at recombination 



net* 



r s (a*) ' 

where r(z) = j dx j H(x) is the comoving distance and r s (a*) is the comoving sound 
horizon at recombination 

pa * c s (a) 



a 2 H{a) 



da, 



where the sound speed c s (a) = l/y3(l + i?&a) and Rb = 30^/4fi^ is the photon-baryon 
energy density ratio. The x 2 of the CMB data is constructed as: 



where 



X 



XTC M,cmb X ' 



( l a - 302.09 
R- 1.725 
yz* - 1091.3 



(3.6) 



and the inverse covariance matrix 



^M.cmb 



( 2.305 29.698 -1.333 
29.698 6825.270 -113.180 
-1.333 -113.180 3.414 



V 

The fourth set of observational data is 12 Hubble evolution data from [26] and [27]. 
Its Xh ls defined as 



12 



[H(zi) - H ob ( Zi )i< 



i=l 



Note that the redshift of these data falls in the region z E (0, 1.75). 

The Big Bang Nucleosynthesis (BBN) data we use here are from [28, 



2 (n b0 h 2 - 0.022) 



0.002 2 



(3.7) 

whose x 2 is 
(3.8) 
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[40] 


0.77 


0.91 


0.36 


[44] 


0.78 


0.70 


0.08 


[40] 


1.4 


0.90 


0.24 
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3.0 


1.46 


0.29 


[46] 



Table 1: Currently available data for linear growth rate f b s used in our analysis, z is rcdshift; 
a is the la uncertainty of the growth rate data. 



And finally for the growth factor data, we define 



x , = ^M_W!f. (3 . 9) 

i=i 1 

The 11 data of growth factor are summarized in Table 1 [38]. 
3.3 Fitting Results 

The best fitting values and errors of the model parameters are summarized in Table 2, 
where we also list the best fitting values of the corresponding parameters of ACDM model 
for comparison. The best fitting values of f2 m o and h are slightly smaller than corresponding 
ones in the ACDM model and the best fitting values of £lbo are larger than corresponding 
ones in the ACDM model. We also can see from Table 2 that adding the data of growth 
factor dose not have much impact on the values of the parameters, both at la confidence 
level and 2a confidence level, which may mean that this model is not sensitive to the 
linear growth rate of matter. In addition, we find that 7 < 1 is excluded at 2a confidence 
level, which means that in the ghost dark energy model, (3 < 0. Furthermore, we see that 
the subleading term H 2 of the dark energy density, the early dark energy, could have a 
fraction energy density around 10%. In Figure 1 and Figure 2, we plot the ID marginalized 
distribution probability of each parameter using the full datasets. 

In Figure 3 we plot the evolution behaviors of the equation of state w(z) of DE and 
the deceleration parameter q(z), with the best fitting values of the model and the ACDM 
model. In the calculation, we employ the following relations: 

w{z) = - l+ mz)^z-> (3 - 10) 



-7- 



parameter 


SN+BAO+CMB+H+BBN 


SN+BAO+CMB+H+BBN+F 


KCDM 


h 


n Ml'/ / 1 ( i ' i r. 1 

U - D4Z +0.012, +0.023 


n in - . 1 \ / w -' 

n f?/i q — u.uio, — u.uz ; 
U - D4/: +0.010, +0.021 


0.708 




U - ZDU +0.014, +0.026 


„ 9C ... -0.014, -0.025 
U.zol_|_o.oi3, +0.026 


0.266 




u - uoz +0.002, +0.003 


U.UOz^o.002, +0.003 


0.045 


7 


,,,-0.035,-0.062 
i - ii4 +0.029, +0.058 


i - iuo +0.035, +0.063 


\ 


fo 


\ 


n -70-0.018,-0.029 
U - 4M +0.012, +0.024 


0.485 



Table 2: The best fitting values within la and 2a errors for h, fl m o, flbo, 7 and /o for the dark 
energy model. The second column shows the results using the datasets without the data of growth 
factor, and the third column shows the results fitted with full datasets. The last column shows the 
best fitting results of KCDM model using the full datasets for comparison. 




0.60 0.62 0.64 0.66 0.68 0.22 0.24 0.26 0.28 




0.048 0.050 0.052 0.054 0.056 1.05 1.10 1.15 1.20 



Figure 1: ID marginalized distribution probability of h, H m o, flbo and 7 using the full datasets. 



The results show that in the QCD ghost dark energy model, the universe transits from early 
matter dominant phase to the de-Sitter phase in the future, as expected. The accelerating 
expansion begins at z = 0.75, which is earlier than what the ACDM model predicts. w(z) 
varies from w > — 1 to w = — 1 which is similar to freezing quintessence model [47]. 

The total x 2 of the best fitting values of this model using the full datasets is Xmin = 
589.422 for dof = 586. The reduced x 2 equals to 1.006, which is acceptable, but Xmin 1S a 
little larger than the one for the KCDM model, Xacdm = 558.890. A similar conclusion 
is also reached by other authors using different data set in [48]. That work studies the 
dynamics of varying vacuum energy as a cosmological constant. That is, the equation of 
state of the vacuum energy is always kept as w = — 1. In that case, there must exist 
interaction between matter and the vacuum energy. In Figure 4, we plot the evolution of 
the growth factor for the QCD ghost dark energy model and the KCDM model, it shows 
that the ghost dark energy model can not be discriminated by the data, and that both of 
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0.44 0.46 0.48 0.50 



fo 

Figure 2: ID marginalized distribution probability of fo using the full datasets. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



Figure 3: Evolution behaviors of the equation of state of DE and the deceleration parameter for 
the QCD ghost dark energy model and the AC DM model. 

these two models fit the data very well, even the ghost dark energy model looks fitting the 
data better. 

4. Conclusion and Discussion 

The accelerating expansion (dark energy) of the universe must be closely related to the 
vacuum energy of quantum fields. It is believed that the difference between the vacuum 
energies in Minkowski space and in FRW universe might be the origin of observed dark 
energy. However, the naive estimate indicates that the difference should be of the form 
H 2 Ac [20, 21]. Such a term is too small and cannot derive the universe to accelerating 
expansion. But this term may play an important role in the early evolution of the universe, 
acting as an early dark energy. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

z 

Figure 4: Evolution behaviors of the growth factor for the QCD ghost dark energy model and the 
ACDM model. 

On the other hand, the vacuum energy difference from the Veneziano ghost field intro- 
duced in order to solve the so-called U(1)a problem in QCD has the exact form, aH + f3H 2 , 
where a ~ 

A QCD ~ (WOMeV) 3 . The leading term gives exactly the order of the observed 
dark energy. Therefore the QCD ghost dark energy model is very attractive in the sense 
that this model needs not introduce new degrees of freedom or modify Einstein's general 
relativity, to explain the accelerating expansion of the universe observed today. 

In this paper, based on the vacuum energy of QCD ghost field, we investigated a DE 
model whose energy density has the form aH+f3H 2 . We studied the dynamical evolution of 
the QCD ghost dark energy model and fitted this model with observational data including 
Snla, BAO, CMB, BBN, Hubble parameter and the growth factor. The best fitting results 
show that the subleading term of the energy density makes a negative contribution to the 
total energy density. In this model, the universe transits from early matter dominant phase 
to a de-Sitter phase in the future, and the accelerating expansion begins at z = 0.75, which 
is earlier than that of ACDM model. The equation of state of DE varies from w > —1 to 
w = —1 like a freezing quintessence model. 

The total x 2 OI the best fitting values of this model is x m in = 589.422 for the full 
datasets with dof = 586. The reduced x 2 is 1.006, which is acceptable, but x m in 1S a little 
larger than the one for the ACDM model, Xacdm = 558.890, for the same datasets. We 
further studied the cosmological dynamics of the model by considering the effect on the 
growth rate of matter. The ghost dark energy model can not be discriminated by the data, 
and both of this model and the ACDM model fit the data very well. 

Finally before ending this paper, we would like to stress that in fact there have not 
been any precise calculations showing that the vacuum energy density of the Veneziano 
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ghost of QCD in a FRW universe is of the form, aH + /3 2 , because the vacuum energy 
calculation of the Veneziano ghost is quite difficult in both flat and curved spacetimes due 
to the intrinsic difficulties of QCD and strongly interacting fields in general [10]. But the 
vacuum energy calculations of the Kogut-Susskind ghost in 2d QED (which is the direct 
analogue of the Veneziano ghost in QCD), in 2d topological nontrivial spacetime and curved 
space [10, 16], and the vacuum energy calculations of the effective Veneziano ghost of QCD 
in 4d Rindler spacetime [11] indeed indicate such a power-law behavior. Note that in 
those calculations, the kinetic contribution is not included, it is expected that the kinetic 
contribution is of the same order of magnitude of the potential. While the dynamics of the 
effective scalar field model for the Veneziano ghost of QCD in FRW universe is discussed 
in the fourth reference in [10], here we have further made an approximation that the size 
scale L of some nontrivial manifold is replaced by the Hubble size of the FRW universe, 
which becomes time-dependent. In principle, in such an extension, one has to consider the 
kinetic contribution of the time-dependent scale to the vacuum energy density. However, 
we assume that such a contribution is subdominant due to slowly evolution of the scale, 
compared to the potential term, and absorb some uncertainties into the two coefficients 
a and f3, because at the moment one could not make an explicit calculation to take into 
account the effect. Of course, the kinetic contribution of the Veneziano ghost of QCD and 
the effect in a time-dependent spacetime should be seriously studied if the QCD ghost dark 
energy model studied in this paper can fit well with the observational data. Clearly at the 
moment one just can review this model at a phenomenological level. 
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